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We analyze two pairs of doublet states for the two Anderson impurity problem. We found that 
for short interimpurity distances they have a lower energy than the ferro triplet and the antiferro 
singlet. For medium to long distances between the impurities, the doublets also have a lower energy 
than the two decoupled Kondo-singlet state for most of the range of interest. The mechanism behind 
their energy gain is an coherence enhanced Kondo-like interaction. 
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The two Anderson impurity problem, i.e. how two 
magnetic impurities interact when they are embedded in 
a metallic host, has been the subject of many theoret- 
ical studies, ranging from perturbation theoryi^^ and 
narrow band approximation^ to renormalization group 
analysis 5 :. Besides its natural presence in alloys, today 
this system can be tailor made by means of quantum 
dots&^ and promises to be a relevant circuit component 
in quantum electronics. Here we present a variational 
wave function analysis that, as in the single impurity 
case&, gives a clear description of the system. When 
needed, we specialize in the one dimensional (ID) case 
because of its technological applications. 

The Anderson Hamiltonian for magnetic impurities di- 
luted in a metallic host is: 

H = E ^aCka + ^(e*^ d) a C ka + h.C.) 
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where the fermion operators c ka act on the conduction 
band states and dj a on the orbital of the magnetic im- 
purity situated at rj. Single state energies e k ,—Ed are 
refereed to the Fermi energy (Ep — hp /2m), i.e. there 
is an implicit —fjt N term in the Hamiltonian that regu- 
lates the population of the system (/x = Ep is the chem- 
ical potential and N the total number operator), V is 
the d — c hybridization and N c is the number of cells in 
the metal. In the Kondo limit the impurity levels are 
well below the Fermi energy (Ed 3> 0), and they can 
not be double occupied due to the Coulomb repulsion in 
them (U 3> 2E d ). In order to simplify calculations we 
renormalize the vacuum (denoted by \F)) to be the con- 
duction band filled up to the Fermi energy and we make 
an electron-hole transformation for band states below the 
Fermi level: b\ G = c_fc._ CT for \k\ < hp . In this way the 
energy of a hole excitation (i.e., to remove an electron 
from below the Fermi level) is explicitly positive. 

We consider here the two impurities case, one impurity 
at — r/2 (the Left impurity) and the other at r/2 (the 
Right impurity). Although calculations were made with 
the second quantization operators, we use in the text a 



"ket" notation for the impurity states: the first symbol 
indicates the population status of the Left impurity and 
the second the status of the Right impurity, e.g. | J,, |) = 
d[/ m \F), |0,T) = 4 T |F), |0,0) = |F),etc. 

In this situation (Ed 3> 0, U 3> 2E d ) the lower energy 
configurations of the system are the ferro triplet (| f, f 
), (I L4) + I 4,1"}), I 44)) and the antiferro singlet 
(I Ti-L) — I 4>T))i W1 th energies —2E4. These states are 
not eigenstates of the system; the hybridization mixes 
these states with configurations having excitations in the 
band. Perturbation theory allows for the calculation of 
the corresponding corrections. These energy corrections 
are the same for both the ferro (FF) and antiferro (AF) 
states (and twice the single impurity corrections) up to 
the well known RKKY-like four order term, which has 
opposite contributions on themA£: 
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where the q (k) sum is over hole (electron) excitations 
(\q\ < k F , \k\ > k F ) and v = V/^fW c . 

A related state is the Kondo singlet, the ground state 
when just one impurity is considered. If only the Left 
impurity is present, this singlet is described by the vari- 
ational wave function^ 
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where the variational parameter results to be — 
v/(Es + Ed — e q ) and a self-consistent equation is ob- 
tained for the singlet energy Es, 
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after the pole in the sum, Eg = — E d — 8k is proposed 
(Ed ^> 5k), and Eq.© becomes an equation for the 
Kondo energy 5k, 
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from which 
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is obtained. Jn (= (V 2 / Ed) n , n Q being the density 
of states at the Fermi level) is the relevant parame- 
ter for these theories, and D, the half band-width, be- 
comes the energy scale {Sk, Si? are proportional to it). 
Jn = 0.07238 gives Sk one thousandth of D, a reasonable 
value. 

We study two doublet states, the even doublet, 

l^) = (IT,o) + |o,t)) + V{ 
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and the odd doublet, 

|£>o) = (|T,0)-|0 ) T)) 
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After a standard minimization procedure, the varia- 
tional parameters are obtained 



E x + 2E d ± E, 
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where Ex stands for the energy of the doublet (Ee or 
Eo) and the plus sign in front of E# holds for the ferro 
a factor (a F ) and the minus for the antiferro. For the 
energies, self-consistent equations result 
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For Eo a similar equation holds, with the sin and cos 
interchanged. This behavior, as well as the sign of E#, 
is due to the spatial symmetry of the ferro and antiferro 
states. As is the case of the Kondo singlet, the dominant 
pole in the sums of Eq. (|llf) determines the value of Ee ■ 

First, let us analyze the situation when E/j vanishes. 
This situation corresponds to the zeros of the RKKY 
interaction and the fcpr 3> 7r limit. 

For En = 0, after expanding the sin 2 , cos 2 functions, 
Eq.lJTIJ reads 
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Proposing Ee = —2E d — Se we obtain 
_ 9 2 — cos q.r 
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Therefore for very large r, such that the contribution 
from the cos term vanishes, the energy gain of the doublet 
(8e) tends to that of one Kondo singlet. This is not 
surprising, given that the doublets, when E# ■ — 0, can be 
written as 



\Dx) = \1l)®\Sr)±\S l )®\U), 
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i.e., the combination of a Kondo singlet in one impurity 
and the other single occupied plus (minus) their mirror 
image. 

This energy, Se, must be compared with twice the 
Kondo energy, because, for very large r, the latter is 
the energy gain for a state that has simultaneously both 
impurities forming a Kondo singlet (\Sl) ® l<Sii))- 

Now, we examine the effect of the coherence term in 
Ea. p3|l : here we present the ID analysis. The sums in 
Eci. (|13|) can be done exactly. The first one is the Kondo 
integral and we call the second one the Doublet Coher- 
ence integral, 

W c ^S~hr = n olK(S) = n ln(l + ^) , (15) 
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where x = kEr. The Coherence integral has a logarithmic 
dependence on S that goes like Ik , so it is useful to define 
Cd(S,x) — Ix(S,x)/Ik(S), which is a decaying oscilla- 
tory function that depends lightly on S. We will see that 
Cd plays the role of a coherence driven extra connectivity 
between the components of the doublet wave function. In 
Fig^we plot Co as a function of kEr for various values 
of S, and the r dependence of the 1D-RKKY is also plot- 
ted. It can be seen that the RKKY reaches its asymp- 
totic form ( — cos (2x)/ttx ) very quickly, at fc^r > n, 
whereas Ix approaches a similar form ( sin (a;) D/xS ) 
but at kEr 8/D > n, a distance about one thousand 
times greater, the Kondo length . 

As a function of Ik and Cd Eq. I|13|) becomes 
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from which Se can be obtained numerically. It is inter- 
esting to write it the Kondo way (Eq.Q), 



Se = D exp 
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which can be taken as an approximate solution by evalu- 
ating Cd at Sk, or used iteratively to find the exact result 
with fast convergency. Note that Cd couples to the "con- 
nectivity" factor (the 2 in the exponent of Eq.Q), and 
small changes in this factor produce huge changes in Se- 
For So the sign in front of Cd is positive. These energies 
are plotted in Fig^Jas a function of k F r for Sk = 0.001 D. 
They alternatively surpass 2 Sk up to k F r ~ 50-7T. 

This effect, Se,So > %Sk, is present for values of Jn 
lower than Jn c — 0.240 (which gives Sk — D/8); the 
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FIG. 1: Range functions for the RKKY, Jn(2kFr), and the 
Doublet Coherence factor, Co(S,kFr). The period of the 
RKKY is half the one of Cd and its amplitude decays very 
rapidly. These characteristics are due to the fact that the 
RKKY depends on two excitations, whereas the Doublet Co- 
herence involves just one. 
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FIG. 2: Correlation energy gains of the odd, So, and even, Se, 
doublets as a function of the distance between the impurities, 
for E_r = and Jn = 0.07238. The inset is Se up to very long 
distances, and it periodically beats 28k up to ~ 507T. 



lower Jn, the greater this effect. The periodicity with 
which \Do), \De) alternate is approximately twice the 
RKKY period (they are ruled by Cd), and they stand 
in this situation up to very long distances. These prop- 
erties are due to the fact that the Doublet Coherence 
integral involves just one hole excitation. At the inter- 
section points one has 5e,So = $K an d Cd = . 

The situation for Ejj ^ can be analyzed with the 
same concepts already introduced. We analyze here the 
even doublet energy Ee when £7? > 0. The analysis for 
Eo and/or Tir < is very similar. 

For Y.r > the lower energy configurations are the 
ones of the ferro triplet, which makes the pole of 
the relevant one in the equation for Ee, Ea. (|ll|) . Thus, 



proposing Ee = — 2£^ — £r— Se, Ea. (|llf) is transformed 
to 
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and then, using the definitions for Ik and Ix, we have, 
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a self-consistent equation for Se (in Ea. (|20(l Ix is also 
a function of kpr, but only the "energy" dependence is 
quoted). If Ed > "Er,5e the left side of Eq.(j20J is just 
1/ Jn; if not, a new parameter must be introduced: the 
Ed/D ratio. In the present work we keep it simple (just 
Jn). For So (and > 0) the sign of the Ix terms are 
inverted. Therefore when So decreases Se increases and 
viceversa, following the oscillations of Cd as a function 
of kpr, just as in the Ej? = case. For < the 
dominant pole is that of the antiferro component of the 
doublets and Se, So swap equations, with |£r| instead of 
E/j and one important difference: the bare connectivity 
of the antiferro component is 1/2, instead of the 3/2 of 
the ferro one, and hence these energy gains are lower than 
in the E# > case. 
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FIG. 3: Total energy gain of the even doublet, Ab = + 
Se, as a function of r. Note that in the relevant ranges for 
this state the main contribution to Ab comes from Se ■ 

In Fig|3| we show Se, ^r, and the total energy gain 
of the even doublet, A# = |Sjj| + Se- Note the kinks 
in Se at the zeros of E_r, those kinks reflect the change 
of dominant component in \De) when T,r changes sign. 
Nevertheless, the total energy gain is smooth at those 
points. In the regions in which Se is equal to zero the 
doublet "collapses": i.e. when Se — the dominant 



diverges for \q\ = kp, fEas. (|8HU|) . ek F — 0) and 
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thus this is the only surviving configuration of the ones 
that compose the doublet, and a FF (AF) state with 
just one hole at the Fermi level is simply a FF (AF) 
state. The kinks of in these regions are because A# = 
|£r|. In the relevant regions (S F > So) usually 5e > 
thus the dominant pole configurations (FF (AF) 
based) have just a little more weight in the doublet than 
the "dominated" one (AF (FF) based), and the = 
regime is soon reached. 

Finally, in FigQ] we plot both and Aq as function 
of r, is also shown (8 k = 0.001Z?). For this value of 
Jn, in the measure that r increases the system alternates 
between the odd and even doublets, just as in the Y,r = 
case, except for a small region above k F r = tt/2. At 
this point is near an AF maximum and Cd is at its 
first zero. Hence there is not connectivity enhancement, 
neither for the odd, nor for the even doublet, and they 
collapse into the AF state. 
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FIG. 4: A_b and A as a function of r, for 5 K = 0.00173. 
In conclusion, we give a clear description of the behav- 



ior of these doublet states of the Two Anderson Impurity 
System and their interplay with the RKKY interaction. 
A remarkable point of the model is the long distances at 
which the impurities remain in a correlated state. A full 
understanding of this property, which our model provides 
in detail, allows these systems to be used in the design 
of quantum devices. 
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Appendix: In ID, S R , Eq.J^Jl, half the RKKY interac- 
tion, is 
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= (- Jn) 2 D J R (2k F r) , (A.21) 
where the r dependence is given byS, 



J R (2k F r) = - [J - Si(2k F r)} . 
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The ID Doublet Coherece integral is given by 



(A.22) 



Ix (6, x) — cos (x + xd) [Gi(x + xd) — Ci(x)] + 

sin (a; + xd) [Si(x + xd) - Si (a;)] , (A.23) 

where d = 8/D. Ci (Si) is the Coslntegral (Sinlntegral) 
function, as defined in Mathematica®. 

The upper critical value of Jn, given that \Cd\ < 1, is 
defined by the equation 



e -l/(3Jn c ) = 2 e -l/(2/n c ) 

which gives Jn c = 1/6 In 2. 
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